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1 Introduction 



In integrable models of statistical mechanics an important object is the 7^-matrix 
TZ{u), where u is the spectral parameter. It acts on the tensor product ® for a 
given vector space V and it is the solution of the Yang-Baxter (YB) equation 

7^12(M)7^13(M + v)n23{v) = n^Mn^ziu + v)ni2{u) (i.i) 

in ^V"^ ® V"^, where IZ12 = 1Z®X^ 7^.23 = 1 ®TZ, etc. and X is the identity matrix. 
If 71 depends on a Planck-type parameter r] so that 7l{u, rj) = l + 2ri r{u) + 0(77^), then 
the "classical r-matrix" obeys the classical YB equation 

[ri2(M), ri3(M + v)+ r23{v)] + [ris^u + v),r23{v)] = 0. (1.2) 

Nondegenerate solutions of (|1.2|) in the tensor product of two copies of simple Lie 
algebra g , rij{u) G gi ® gj , i,j = 1, 2, 3, were classified by Belavin and Drinfeld 0. 

The classical YB equation interplays with conformal field theory in the following 
way: In the skew-symmetric case rji{—u) +rij{u) = 0, it is the compatibility condition 
for the system of linear differential equations 

^ d'^izi zn) ^ Y^rij{zi - Zj)'^{zi,...,ZN) (1-3) 

in N complex variables zi, ...,zn for vector- valued functions with values in the tensor 
space V = ^ ■ ■ ■ ® V'^ and k is a coupling constant. 

In the rational case 0, very simple skew-symmetric solutions are known: r{u) = 
C2/U, where C2 € g ® g is a symmetric invariant tensor of a finite dimensional Lie 
algebra g acting on a representation space V. Then the corresponding system of lin- 
ear differential equations (|1.3| ) is the Knizhnik-Zamolodchikov (KZ) equation for the 
conformal blocks of the Wess-Zumino-Novikov-Witten (WZNW) model of conformal 
theory on the sphere 0. 

The algebraic Bethe Ansatz is the powerful method in the analysis of integrable 
models. Besides describing the spectra of quantum integrable systems, the Bethe 
Ansatz also is used to construct exact and manageable expressions for correlation 
functions p. Various representations of correlators in these models were found by 
Korepin using this method. 

Recently, Babujian and Flume developed a method which reveals a link to the 
algebraic Bethe Ansatz for the theory of the Gaudin model. In their method the wave 
vectors of the Bethe Ansatz equation for inhomogeneous lattice model render in the 
semi-classical limit solutions of the KZ equation for the case of simple Lie algebras of 
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higher rank. More precisely, the algebraic quantum inverse scattering method permits 
us write the following equation 

t{u\z)^{ui^...Up) = A(m, Ml, • • • , Up\z)^{ui, ■■■,Up) -V] . (1.4) 

Here t{u\z) denotes the transfer matrix of the rational vertex model in an inhomo- 
geneous lattice acting on an iV-fold tensor product of SU{2) representation spaces. 

= ■ ■ ■ Ua-i,u, Ua+i, Up). J-'a{ui, ■ ■ ■ , Up\z) and A{u, Ml, ■ ■ ■ , Up\z) are c num- 
bers. The vanishing of the so-called unwanted terms, J-'a = 0, is enforced in the usual 
procedure of the algebraic Bethe Ansatz by choosing the parameters ui, ...,Up. In this 
case the wave vector $(mi, ■ ■ ■ , Up) becomes an eigenvector of the transfer matrix with 
eigenvalue A(u, ui, ■ ■ ■ , Up\z). If we keep all unwanted terms, i.e. J-'a 7^ 0, then the wave 



vector $ in general satisfies the equation (1^), named in as off-shell Bethe Ansatz 



equation (OSBAE). There is a neat relationship between the wave vector satisfying the 
OSBAE (|1.4| ) and the vector- valued solutions of the KZ equation (|1.3|) : The general 
vector valued solution of the KZ equation for an arbitrary simple Lie algebra was found 
by Schechtman and Varchenko . It can be represented as a multiple contour integral 

\l/(zi, . . . , zn) ^ " ' ^{ui, Up\z)(f){ui, ...,Up\z)dui ■ ■ ■ dup. (1.5) 

The complex variables zi,...,zn of (|1.5|) are related with the disorder parameters of 
the OSBAE . The vector valued function 0(mi, ...,Up\z) is the semi-classical limit of the 
wave vector $(ui, Up|z). In fact, it is the Bethe wave vector for Gaudin magnets |jTO|, 
but off mass shell. The scalar function X{ui., ...,Up\z) is constructed from the semi- 
classical limit of the K{u = z^^ui, ...,Up\z) and Taiui, ■ ■ ■ ,Up\z). This representation of 
the A^-point correlation function shows a deep connection between the inhomogeneous 
vertex models and the WZNW theory. 

In this work we generalize previous results applying the Babujian-Flume ideas for 
osp{2\l) rational solution of the graded version of the YB equation [llj. It is shown 



that this ideas persists for the case of semi-classical limit which corresponds to the 
osp{2\l) trigonometric r-matrix. 

The paper is organized as follows. In Section 2 we present the algebraic structure 
of the osp(2|l) vertex model. Here the inhomogeneous Bethe Ansatz is read from the 
homogeneous case previously known |l^. We also derive the Off-shell Bethe Ansatz 
equation for the fundamental representation of the osp(2|l) algebra. In Section 3 , 
taking into account the semi-classical limit of the results presented in the Section 2, we 
describe the algebraic structure of the corresponding Gaudin model. In Section 4, data 
of the off-shell Gaudin equation are used to construct solutions of the trigonometric 
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KZ equation. In Section 5, our results are extended for the highest representations of 
the algebra osp(2|l). Conclusions are reserved for Section 6. 

2 Structure of the osp(2|l) Vertex Model 

We recall that the osp(2|l) algebra is the simplest superalgebra and it can be viewed 
as the graded version of s/2. It has three even (bosonic) generators H, X"^ generating 
a Lie subalgcbra s/2 and two odd (fermionic) generators , whose non-vanishing 
commutation relations in the Cartan-Weyl basis reads as 

[H,X^] = ±X±, [X+,X-]^2H, 

[H,V^] = ±^V^, [X^,V^]^V^, [X^,V^]^0, 

{V^,V^} = ±^X±, {V\V-}^~H. (2.1) 
The quadratic Casimir operator is 

C2 = H' + X-} + [V+, V-], (2.2) 

where {• , •} denotes the anticommutator and [• , •] the commutator. 

The irreducible finite-dimensional representations pj with the highest weight vector 
are parametrized by half-integer s — j/2 or by the integer j G A^. Their dimension 
is dim(pj) = 2j -|- 1 and the corresponding value of C2 is j{j + l)/4 — s{s + 1/2), 
8 = 0,1/2,1,3/2,... 

The representation corresponding to s = is the trivial one-dimensional represen- 
tation. The s > 1/2 representation contains two isospin multiplets which belong to 
isospin s and s — 1/2, denoted by \s,s,m) and \s,s — 1/2, m), respectively. The first 
quantum number characterizes the representation and the second and third quantum 
numbers give the isospin and its third component. After a convenient normalization 
of the states , a given s-representation is defined by 

H\s,s,m) = m\s, s,m) , 
X'^lSjSj'm) — ^ {s ^ m){s ± m + 1) \s, s, m ± 1) , 

±-^{s=fm) \s,s- 1/2, m± 1/2) , 

H\s,s — 1/2, m) — m |s, s — 1/2, m) , 
X^|s,s-l/2,m) = ^/(s^^T72^m)(7^T72T7^H^|s,s-l/2,m±l), 

l/2,m) = ±^\/(s - 1/2 ± m + 1) |s, s, m ± 1/2) . (2.3) 
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The fundamental representation has s = 1/2 and is given by 



H 





(2.4) 
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In (12^41) the basis is 
considered as even and the second as odd, i.e., our grading is BFB. 
In the j-representation the odd part has the form |]l3l : 



\-,\-,—\)- The first and third vectors will be 
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(2.5) 



where 

(r,_i,\/,_2,^,-3,...,V^-,) 

(iy„iy,_i,PF,_2,...,iy-,+i) 

For the even part we can see from 
due to isospin integer: 



1 



l,v^,...,v^, 



(- Jj, yi, -J J - 1, J2,..., -v^, Jj). (2.6) 



that H is diagonal and always has eigenvalue 



H = -diag(j,j - 1, 1,0, -1, -j). 



(2.7) 



Moreover, are given by the s/2 composition which results in a clear relation with 
the odd part: X± = ±A{V^f. 

2.1 OfF-Shell Bethe Ansatz Equation 

Consider V = Vo(BVi a, Z2-graded vector space where and 1 denote the even and odd 

s 

parts respectively. The components of a linear operator A ® B in the graded tensor 

s 

product space V ® V result in matrix elements of the form 



{A ® Byjp = (_)p(/5){f(")+p(7)) A^^Bps 



(2.8) 
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and the action of the permutation operator V on the vector \a) 
given by 



eV®Vis 



V \a) 



\p{a)p{f3) 



I") =^ = i-r^'-^'^'^S^, 5,s, (2.9) 



where p{a) = 1 (0) if \a) is an odd (even) element. 

Besides TZ we have to consider matrices R = VTZ which satisfy 

Ri2{u)R23{u + v)Ri2{v) = R23{v)Ru{u + v)R23{u). (2.10) 

The regular solution of the graded YB equation for the fundamental representation 



of osp{2\l) algebra was found by Bazhanov and Shadrikov in |jT^. It has the form 



R{u,f]) 
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(2.11) 



where 



Xi 

X5 
xe 

1/7 



sinh(M + 2ri) sinh(M + Si]), X2 = sinhMsinh(M + St]), 
sinhMsinh(M + 77), X4{u) = sinhMsinh('U + Sr^) — sinh2r7sinh3r7, 
e~"sinh2r7sinh(M + Srj), = e"sinh2?7sinh(M + Srj), 

— ee~"~^'' sinh 2r] sinh u, ye = ee"'^^'' sinh 2r] sinh u, 
e~"sinh2?7 (^sinh('U + 3?]) + e"'' sinh-ul 
e" sinh 2// (sinh(M + Si]) + e'' sinh u) . 



(2.12) 



where e = ±1. Here we have assumed that the grading of threefold space is p(l) = 
p(3) = and p{2) = 1 and we will choose the solution with e = 1. 

Let us consider the inhomogeneous vertex model, where to each vertex we associate 
two parameters: the global spectral parameter u and the disorder parameter z. In this 
case, the vertex weight matrix TZ depends on u — z and consequently the monodromy 
matrix will be a function of the disorder parameters Zi. 
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The graded quantum inverse scattering method is characterized by the monodromy 
matrix T{u\z) satisfying the equation 



R{u 



T{u\z) ® T{v\z) 



T{v\z) ® T{u\z) 



R{u 



(2.13) 



whose consistency is guaranteed by the graded version of the YB equation (|1.1|). T{u\z) 
is a matrix in the space V (the auxihary space) whose matrix elements are operators on 
the states of the quantum system (the quantum space, which will also be the space V). 
The monodromy operator T{u\z) is defined as an ordered product of local operators 
Cn (Lax operator), on all sites of the lattice: 



T{u\z) = Cn{u - zm)Cn~i{u - zm-i) ■ ■ ■ Ci{u - zi). 



(2.14) 



The Lax operator on the n*^ quantum space is given the normalized graded permutation 
of (iH: 



r = 

'~'n 
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(2.15) 



Note that L^^piu), 
It means that the monodromy matrix has the form 



a, /5 = 1, 2, 3 are 3 by 3 matrices acting on the n*'^ site of the lattice. 



(2.16) 
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where 



Tij{u\z) 



ki, 



1,2,3. 
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E 

■,fc]V-l = l 



l\Z\u - Zn) ® L)!' ''{u - zn-i) 



{N-l). 



L 



(2.17) 
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The vector in the quantum space of the monodromy matrix T{u\z) that is annihi- 
lated by the operators Tij{u\z), i > j {Ci{u\z) operators, i = 1,2,3) and it is also an 
eigenvector for the operators Tii{u\z) ( Ai{u\z) operators, i = 1, 2, 3) is called a highest 
vector of the monodromy matrix T{u\z). 

The transfer matrix t{u\z) of the corresponding integrable spin model is given by 
the supertrace of the monodromy matrix in the space V 

3 

Tiu\z) = Tuiu\z) = Aiiu\z) - A2iu\z) + Asiu\z). (2.18) 

i=l 

Algebraic Bethe Ansatz solution for the inhomogeneous osp{2\l) vertex model can 



be obtained from the homogeneous case [T^. The only modification is a local shift of 
the spectral parameter u u — Zi. 

Here we will define some functions that will be used in the calculations of the Bethe 
Ansatz: 



z{u) 
u{u) 
Ziuk - Uj) 



Xi{u) sinh('u + 2?7) 



X2(U] 



sinh-u 

xi{u)x^{u) 



y{u) 



X'iiu) sinh(M + ?7) 

?/6 {u) QU+2r) ^lYih 2?7 ' 

sinh(u + rf) 



X4{u)xs{u) — XQ{u)yQ{u) smh.{u — 7]) 

z{uk — Uj) if k > j 

z{uk — Uj)uj{uj — Uk) if k < j ' 



(2.19) 



We start defining the highest vector of the monodromy matrix T{u\z) in a lattice of 
sites as the even (bosonic) completely unoccupied state 



|0) 



/ 1 


Vo 



(2.20) 



Using ( p.lTp we can compute the normalized action of the monodromy matrix entries 
on this state 



A,{u\z)\0) = Xi{u\z)\0), aiu\z) \0) = 0, Bi{u\z)\0)^{0,\0)} 
X,(u\z) = 



1,2,3. 



(2.21) 



:=lX2[U-Za) 

The Bethe vectors are defined as normal ordered states \1'„(mi, • • • ,m„) which can 



be written with aid of a recurrence formula 15 



^„(Mi, ...,Un\z) = Si(mi|z)*„_i(m2, ...,M„|^) 
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n n 

-52(mi|z)^— ^^^i^ n 2K-Uj)^„_2(M2,...,Mi,...,nJz) , (2.22) 

j=2 ^1^1 ^i) k=2,k^j 

with the initial condition \I/o = |0) , = Bi{ui\z) |0). Here Uj denotes that 

the rapidity Uj is absent: "^{Uj \z) = \l/(-ui, ...,Uj-i, Uj+i, ■ ■ ■ ,Un\z). 

The action of the transfer matrix t{u\z) on these Bethe vectors gives us the following 
off-shell Bethe Ansatz equation for the osp{2\l) vertex model 

n n j—1 

t{u\z)^^{u,,...,u^\z)=AA{uu...,u^\z)-T.^j^^^^^ 

j=l j=2 1=1 

(2.23) 

Let us now describe each term which appear in the right hand side of (|2.23| ) ( for 
more details the reader can see [|12]): In the first term the Bethe vectors (|2.22| ) are 
multiplied by c-numbers A„ = A„(m,mi, ...,m„|z) given by 

An = X,iu\z) n ziu, -u)~ {-TX2{u\z) \{ '^"^ ~ + XMA n ^4^i^. 

(2.24) 

The second term is a sum of new vectors 



multiplied by scalar functions jFj""^-* given by 

n n 

Tf-^^ = X^{u,\z) Yl Z{uk - u,) + {~rX2iu,\z) l[ Z{uj - Uk). (2.26) 

Finally, the last term is a coupled sum of a third type of vector-valued functions 

*;i-2) = B2{u\z)^n^2{uu «,), (2.27) 

with intricate coefficients 

n 

= GijXi{ui\z)Xi{uj\z) W Z{uk- ui)Z{uk- Uj) 

k=l,k^j,l 

n 

-{-)'^YijXi{ui\z)X2{uj\z) Jl Z{uk- ui)Z{uj - Uk) 

k=i,k=/=j,i 

n 

-{-)'' FijXi{uj\z)X2{ui\z) Jl Z{ui-Uk)Z{uk-Uj) 

k=i,k=/=j,i 

n 

+HijX2{ui\z)X2{u,\z) n Z{u,-Uk)Z{ui-Uk). (2.28) 
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where Gij , Yij , Fij and Hij are additional ratio functions defined by 

xr{ui-u) 1 z{ui - u) xr,{uj - u) 1 

(j^Zi = — -, T—, r + 



H, 



X3{ui - u) yiui - Uj) uj{ui - u) X2{uj - u) y{u - ui) ' 

y7{u - ui) 1 _ y^{u - ui) 1 
X2.{u - ui) y{ui - Uj) X3{u - ui) y{u - Uj) ' 

y^{u-Uj) y5{u - ui) y5{ui - Uj) 



Yij = —, r {z{u-uij 

y[u-Ui) y X2[U-Uj) X2{U - Ui) X2[Ui - Uj) 

y^{u-ui) \yr^{ui-Uj) 1 z{u - ui) 1 

i'li = 7 T \ 7 T—, T + 



(2.29) 



X2{u-ui) \x2{ui - Uj) y{u - ui) uj{u - ui) y{u - Uj) 
y^ju-ui) 1 \ 
X2{u - ui) y{ui - Uj) J ■ 

In the usual Bethe Ansatz method, the next step consist in impose the vanishing of 
the so-called unwanted terms of ( |2.2^ ) in order to get an eigenvalue problem for the 
transfer matrix: 

We impose JFj""^-* = and J^/""^"* = into ( |2.23| ) to recover the eigenvalue problem. 
This means that \E'„(mi, is an eigenstate of t{u\z) with eigenvalue A„ , provided 

the rapidities Uj are solutions of the inhomogeneous Bethe Ansatz equations 

li z{u,-z.) = (-r^ n "-^r^^Mu,-u,), 

a=l k=l, k^j ^y'^k Uj) 

J = l,2,...,n. (2.30) 



3 Structure of the osp(2|l) Gaudin Model 

In this section we will consider the theory of the Gaudin model. To do this we need to 
calculate the semi-classical limit of the results presented in the previous section. 

In order to expand the matrix elements of T{u\z), up to an appropriate order in 7], 
we will start by expanding the Lax operator entries defined in ( |2.15| ): 



4? = J. + 2r^coth(n - z^) + 2r^' (nl + ^ . ] + o(r^^ 

' 2 smh(M — -Zri 



4? = X„-2r^^ ^^^" ^\ +o{n'). 

smh(u — Zn) 

lS' = I„ - 2,coth. «„ + 2,f (nl + f-ll^) + o(„3). (3.1) 
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and for the elements out of the diagonal we have 



-U+Zn 



sinh(M — Zn) sinh(M — Zn) 



-U+Zn 

sinh(M — ' ' ' ^' ■^^ ' sinh(M — z„ 



p U-\-Zn 

4? = 2r/-^7 ,Xn+o{r^% 4^ = 2r/ — -^-^ + o(r^^).(3.2) 



sinh(-u — Zn) " ' sinh(-u — z 

where V± = 2^^, = 2X± and Tt: = 2iy. 

Substituting and (13) into (|2T^ ) we will get the semi-classical expansion for 
the monodromy matrix entries: 

^ r 

y4i(u|2;) = X + 2?7 ^ COth(M — ZajTi-a + ^"ff < X! COth(M — Za) COth(M — Zh)'Ha ® Ti-b 
a=l (a<b 

:tZa Z^ 




-2b 



N 



sinh(M — Za) sinh(M — Zb) 



4e 



2 sinh(M - Za)"^ 

(3.4) 



y43('u|2;) = X — 2?7 coth(-u — Za)'Ha + 47]^ < coth(n — Zg) coth(n — Zb)'Ha g) Hb 

a=l \^a<b 

^ smh(M - Za) smh(M - Zb) \ 
For off-diagonal elements we only need to expand them up to the first order in rj 

N Q~U+Za 

BAu\..) = -B,(.„|.) = -2,,i;-S^i:— yK-+o(,^ 

N pU-Za 

C^{u\z) = C,{u\z) = 2r^Y.-u S^a+o{v'), 

smh(M - Za) 
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N pU-Za 

C,{u\z) = 2vY.-u S^a+^iv')- (3.6) 

^1 sinh(M - Za) 



Therefore, the semi-classical expansion of the transfer matrix ( p.l8| ) has the form 



t{u\z) =I + 8T]^T^^\u\z)+o{r]^), (3.7) 



where 



^ ' ^ .tle«-sinh(«-.J+„^A2 ^^^4 



sinh(M — Za)2 



(3. 



with 



If s 

Ga{u) = V — — <^ COsh(M - Za) COsh(u - Zb)Ha ® Hb 

^^^smh{za- Zb) I 

+ (^e-""+"''V+ ® - e''^~"'V~ ® Vb+) | . (3.9) 
Here we have used the symmetry 

Qbaiu, Za, Zb) = V Sabiu, Za, Zb) V = Gab{u, Zb, Za), (3.10) 



and the identity 

/ Q-U+Za 



sinh(ti — Za) sinh(M — Zb) sinh(2;a — Zb) I sinh(u — Za) sinh(u — Zb) J 



(3.11) 



The Gaudin Hamiltonians are defined as the residue of t{u\z) at the point u = Za- 
This results in non-local Hamiltonians 

N T , - 1 



Ga = E-T7 .\coshiza-Zb)na^nb + -[e-^^^'>'X^®X,- 

^ smh(2;„ - Zb) [ 2 V 

a = 1,2,...,N. (3.12) 



satisfying 

AT 

S^" °' dzb dz, 



Y.Ga = 0, —^ = -^, [Ga,Gb] = 0, \/a,b. (3.13) 



11 



In the next section we will use the data of the algebraic Bethe Ansatz for the 
osp{2\l) vertex model to find the exact spectrum and eigenvectors for each of these 
— 1 independent Hamiltonians. 

Before doing this, we would like to consider the semi-classical limit of the funda- 
mental commutation relation (|2.13| ) in order to get the osp{2\l) Gaudin algebra: 
The semi-classical expansions of T and R can be written in the following form 



T{u\z) = 1 + 2ril{u\z) + o{7f), R{u) =V 1 + 2rir{u) + o{7f 



From (|3 .31 - 13.61) we can see that the "classical /-operator " has the form 

l{u\z) = 



niu\z) -V'{u\z) X-{u\z) 
V+iu\z) V-iu\z) 
X+{u\z) V+iu\z) -n{u\z) 



(3.14) 



(3.15) 



where 



n{u\z) 

V-{u\z) 
X-{u\z) 



N 

J2 COth(u - Za)'Ha, 
a=l 

N ^-U+Za 

-V-. 



N 



E 



^ sinh(u - Za 

N „-u+z„ 



VHu\z) = E 



^ sinh(M - Za 



n 



N 



E 



X^u\z) = E 



^-^smh{u-Za) " 



(3.16) 



coshM 7^ ® 7^ + - (e-" A"- 



e" X 



X- 



(3.17) 



sinh(u - Za 

The classical r-matrix has the form 

r{u) = — ^ 
smhw 

+ e"" {nv^ + v+n) ® {v-n + nv-) 
- e" {V'H + nv-) ® {nv+ + V+H)^ 

Here we observe that we are in the fundamental representation of the osp{2\l) algebra 
where the relation 

{nv^ + v+n){v-n + nv-) = v+v- (3.18) 

holds. Therefore, (|3.17| ) is equivalent to that r-matrix constructed out of the quadratic 
Casimir in a standard way |jT6[. Indeed it corresponds to the second regular solution 
(e = —1) presented in ( |2.12| ). 

Substituting (|3.17|) and ( p.l5|) into the fundamental relation ( p.l3|) , we have 



7^/(1*1^) O l{v\z) + Vr{u - v) 



l{u\z) (g) 1 + 1 (g) l{v\z) 



l{v\z) (g) l{u\z)V + 



l{v\z) + l{u\z) 



Vriu — v) 



(3.19) 
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whose consistence is guaranteed by the graded classical YB equation. 

From ( p. 19 ) we can derive (anti) commutation relations between the matrix elements 
of l{u\z). They are the defining relations of the osp{2\l) Gaudin algebra : 



[n{u\z),n{v\z)] 

[V^{u\z),X^{v\z)] 
[V^{u\z),X^{v\z)] 

[X-{u\z),X+{v\z)] 

{V-{u\z),V+{v\z)} 

[n{u\z),V^{v\z)] 

[n{u\z),X^{v\z)] 

{V^{u\z),V^{v\z)} 



0, 

[X^{u\z),X^{v\z)] = 0, 
2p±(«-f) 

V^{u\z) -V^{v\z) 



sinh('U — 


v) L 






sinh(u — 








sinh(u — 




1 




sinh(-u 


-v) 






sinh(M 


-v) 


1 




sinh(-u 


-v) 



[n{u\z) - n{v\z)] , 
[n{u\z) -n{v\z))], 

e^("~^)V^(M|z) - cosh(M - v)V^{v\z) 
e±(^~-)X^(u\z) - e^^'^'^^X^ivlz))] . (3.20) 



A direct consequence of these relations is the commutativity of t^'^\u\z) 



T^^\u\z),T^'^>{v\z)]=0, Vm,W 



-(2) 



(3.21) 



from which the commutativity of the Gaudin Hamiltonians Ga follows immediately. 
3.1 Off-shell Gaudin Equation 



In order to get semi-classical limit of the OSBAE ( |2.23|) we first consider the semi- 
classical expansions of the Bethe vectors defined in (|2.22| ), (|2.25| ) and (p.27|) : 



^„(ui,...,u„|z) = (-2r7)"$„(Mi,...,u„|2;) +o(r/"+^), 



where 



(n-l) 



2(-277) 



n+1 



V-{u\z)<^n-i{uj |2)+o(r/"+2) 



sinh(M — Uj) 



A A 



fn-2) = -{-^Vr-'X-{u\z)<^n-2{Ul,U,\z)+0{vn 



(3.22) 



$n(wi, •••,Mnk) 



V {Ui\z)^n-liu2, ■■■,Un\z) 

-A'-(^iil^) v4tt^ ^ 

jr^ smh(Mi — Uj] 



$n-2(Mj \Z) 



(3.23) 
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with $0 = |0) and ^i{ui\z) = V"(mi|z)$o- 

Here we would like make a few comments on the structure of these vector-valued 
functions. In (|3.23|) they are written in a normal ordered form. Since we are working 
with fermionic degree of freedom, the function ...,Un\z) is totally antisymmetric. 



$„(Mi, Ui+i, Ui, Un\z) = -$n(^l, Mj, U 



,Un\z). 



(3.24) 



To see this one can use the Gaudin algebra ( |3.2CI| ). For instance, in its antisymmetric 
form the Bethe vector $2 reads as 



$2(^1, M2k) = ^[V {ui\z)V {U2\z)-V {u2\z)V {ui\z)]<!>o 



1 e"i-"2A'-(Mi|z) + e-"i+"2;t'-(M2k) 



l$n. 



(3.25) 



2 sinh(Mi — U2) 

Now we will consider the semi-classical expansions of the c-number functions pre- 



sented in the OSBAE (|2]2|) 



A. 



l + 2(-2r^)2Ai2) + o(r/3), J^l 



(n-l) 



where 



and 



sinh('Uj — ui) \ sinh('U — Ui 



-ft'' 



^U — Ul 



sinh(M — Uj"^ •' 



(3.26) 



N 



N n 



coth(u — Za) coth(M — 



Ml 



a=lj=l 



N 



+ ^ coth(u — Za) coth(u — Zb) + coth(u — Uj) coth(M — Uk), (3.27) 

a<b j<k 



■I 3 



(n-l) 



N 



— y^ COth(Mj — -2a) + X! COth(Mj — Uk) ■ 



(3.28) 



a=l 



Substituting these expressions into the ( |2.23| ) and comparing the coefficients of the 
terms 2(— 2?7)"+^ we get the first non-trivial consequence for the semi-classical limit of 
the OSBAE: 



T^^\u\z) $„(«!,. ..,M„|z) = A(f) $„( 



n[Ui,...,Un\z) . -.f y 



(3.29) 
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Note that in this hmit the contributions from and ^('^_2) are combined to give 

a new vector valued function 

^ n / \k QUj—Uk ^ ^ 

Q\n-l) = V~{U\Z) \z) - X-{U\Z) ■ u/„ _„ X ^n-2{ Uj,Uk \z), 

k=i, k^j ^"i^HUj Uk) 

(3.30) 

where k' = k + 1 for k < j and k' = k for k > j. Therefore, our graded OSBAE 
( ^■29| ) is very similar to that presented by Babujian and Flume for simple Lie algebras 
]T7| . This result could be expected since the superalgebra osp{2\l) has many features 
which make it very close to the Lie algebra [|T^. 

Finally, we take the residue of ( |3.29| ) at the point u = Za to get the off-shell Gaudin 
equation: 

Ga^n{Ul,---,Un\z) = ga^n{Ul, ■ , Un\z) - 



j^J^ giii Za sinh(u; — Za) 

1,2,...,N (3.31) 



where ga is the residue of A^^^ 



N n 

ga = Tesu=za-^n^ = 12 coth(z„ - Zb) -J2 coth(2;„ - ui), (3.32) 

b=/=a 1=1 



and 0(n-i) is the residue of ©(„_!) 



(— ) e^ 

= reS,=..0jn-l) = ^a^n-l(u, \z) - _ . ^n-2{Uk,U, \z). (3.33) 



In this way we are arriving to the main result of this paper. The equation (|3.31 



permits us solve one of the main problem of the Gaudin model, i.e., the determination 
of the eigenvalues and eigenvectors of the commuting Hamiltonians Ga ( p.9| ): ga is 
the eigenvalue of Ga with eigenfunction $„ provided ui are solutions of the following 
equations /j""^'' = 0, i.e.: 

n N 

coth(Mj -Uk) = Y coth(Mj - Za), j = 1, 2, n. (3.34) 

k^j a=l 

Moreover, as we will see in the next section, the off-shell Gaudin equation (|3.31|) pro- 
vides solutions for the differential equations known as KZ equations. 
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4 Knizhnik-Zamolodchickov equation 



The KZ differential equation 



= H,{z)^{z), (4.1) 



appeared first as a holonomic system of differential equations on conformal blocks in 
a WZW model of conformal field theory. Here "^{z) is a function with values in the 
tensor product Vi ® ■ ■ ■ (S> of representations of a simple Lie algebra, k = k + g , 
where k is the central charge of the model, and g is the dual Coxeter number of the 
simple Lie algebra. 

One of the remarkable properties of the KZ system is that the coefficient functions 



Hi{z) commute and that the form uj = J2i Hi{z)dzi is closed |jT9 



In this section we will identify Hi with the osp{2\l) Gaudin Hamiltonians Ga pre- 
sented in the previous section and show that the corresponding differential equations 
( [4.1|) can be solved via the off-shell Bethe Ansatz method. 

Let us now define the vector- valued function \E'(zi, z^) through multiple contour 
integrals of the Bethe vectors ( p.23|) 

^^(^1, ...,Zf^) = j ■ ■ ■ j X{u\z)^n{u\z)dui...dun, (4.3) 

where X {u\z) = X (-ui, ...,Un, zi, zn) is a scalar function which in this stage is still 
undefined. 

We assume that "^{zi, zjq) is a solution of the equations 

«:^^^%f^ = G,vl/(^,,...,^^), a = 1,2,..., AT (4.4) 

where Ga are the Gaudin Hamiltonians ( p.9|) and k is a constant. 
Substituting (O) into (O) we have 



^d^{z^^_^ = / U-^^^M^) + .X{u\zf-^4^] du, (4.5) 



dZa J [ dZa dZa 

where we are using a compact notation § {o} du = §...§ {o} dui ■ ■ ■ dun- 

Using the Gaudin algebra ( |3.20 ) one can derive the following non-trivial identity 

^ . y-,_,. A f °""-^'<»-» ) , (4.6) 

dza dui \sm\i{ui - Za) J ' 
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which allows us write ('J4.5I) in the form 



-ui-\-za ±1 



+- ti-y / ^ (4.7) 
J dui\ smh(M/ - / 

It is evident that the last term of ( |4.7] ) is vanishes, because the contours are closed. 
Moreover, if the scalar function X{u\z) satisfies the following differential equations 



K 



dX{u\z) , , , , dX{u\z) pfn-i'i,,/ , X /,^x 

g^' ' = gaX{u\z), K ' = ff' ^^X{u\z), (4.8) 

we are recovering the off-shell Gaudin equation ( |3.31D from the first term in ( |4.7| ). 

Taking into account the definition of the scalar functions /j"~^'' (pr^ ) and Qa ( ^.32 ), 
we can see that the consistency condition of the system ( [4.8[ ) is insured by the zero 
curvature conditions df'f'~^^ jdza = dga/duj. Moreover, the solution of ( ^4.8| ) is easily 
obtained 

N n N n 

X{u\z) = n sinh(^, - z,f/^ n sinh(M, - Ukf'^ J] 11 sinh(za - u,Y^'\ (4.9) 

a<b j<k a=l j=l 

This function determines the monodromy of "^Izi, z^) as solution of the trigono- 
metric KZ equation ([4.4| ) and these results are in agreement with the Schechtman- 
Varchenko construction for multiple contour integral as solutions of the KZ equation 
in an arbitrary simple Lie algebra 0. 



5 Highest Representations 

The generalization of our results for the highest representations of the osp{2\l) algebra 
requests the knowledge of the corresponding algebraic Bethe Ansatz which. It can 
be obtained from the Bethe Ansatz for the fundamental representation using a fusion 
procedure. Nevertheless, we can use the fact that the one-parameter operator families 
( ^.l(j| ) form the highest weight module of the infinite-dimensional Lie superalgebra. As 
in the SI2 case presented by Sklyanin in [^, it is characterized by the vacuum |0) 

n{u\z)\0) = h{u\z)\0) , V+(m|z) |0) = 0, X+{u\z) \0) = 0, (5.1) 

the dual vacuum (0| 

{0\n{u\z) = h{u\z) {0\, (0| V^(m|z) = 0, {0\ X' {u\z) = 0, (0 | 0) = 1, 

(5.2) 
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and the highest weight scalar function 

N X 3 

h{u\z) = J2'^^ a coth{u- Za), = -, 1, -, ■ • • (5.3) 

a=l 

The solution of the Gaudin eigenvalue problem is given in term of Bethe vectors 
...,Un) given by ( p.23| ) with the operators written in the j-representation ( p.5|) 
and Q. 

It is convenient to express the Gaudin Hamiltonians Ga in terms of a generating 
function t{u\z) which is obtained as the result of replacing H, X"^ and V"^ in the 
quadratic Casimir ( p.2|) by Ti.{u\z), X'^^ulz) and V^(m|z), respectively: 



4t{u\z) = f: — — ^ +j:(ni+ 

t{ e"--^ smh(u - z,) „tl V 



Sa{Sa + 1/2) 



(5.4) 



sinh(M — Za^ ^ 

where Qa{u) is given by ( p.9|) . 

Using a theorem |21|, |22[ which states that $„(«!,...,«„) is an eigenvector of the 
commuting operators Qa{u) or, equivalently, of t{u\z) if and only if the parameters 
Ui, ...,Un satisfy the Bethe equations 

JV n 

2Sa COth('U — Za) = COth('U — Uk). (5.5) 
a=l k=l 

The corresponding eigenvalue A(u) of 4:t{u\z) is then 

A{u) = \^{u)-2du\{u) (5.6) 

where 

N n 

\{u) = ^ 2Sa COth('U — -2a) — X! COth('U — Uk) (5.7) 
a=l fc=l 

From these results we conjecture that the off-shell Gaudin equation for the highest 
representations is also given by ( p.31 ) with 

N n 

cin-1) 



and 



fj = ~ 51 COth(Mj - Za) + Y1 COth{Uj - Uk), (51 

a=l k^j 



ga = Yl 4-5aSfe COth(2:a - Zb) - ^^a COth(2:a - Ul) . (5.9) 
b=/=a 1=1 



Consequently, the monodromy of the function \I^(zi, Z]\j) ( [4.3|) for the j-representation 
is given by 

N n N n 

X{u\z) = n smh{za - Zb)'^'^'"/'' n sinh(Mj- - u^)^/'^ J] 11 sinh(z„ - Uj)''^'^/'', (5.10) 

a<b j<k a=l j=l 

which also is in agreement with the Schechtman-Varchenko construction 
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6 Conclusion 



In this paper a graded 19-vertex model was used to generalise previous rational results 
connecting Gaudin magnet models and semi-classical off-shell Bethe Ansatz of vertex 
models. 

Using the semi-classical limit of the transfer matrix of the vertex model we derive 
the trigonometric osp{2\l) Gaudin Hamiltonians. The reduction of the off-shell Gaudin 
equation to an eigenvalue equation gives us the exact spectra and eigenvectors for 
these Gaudin magnets. Data of the off-shell Gaudin equation were used to show that 
a Jackson-type integral ( ^.3| ) is solution of the trigonometric KZ differential equation. 

In fact, this method had already been used with success to constructing solutions 



of trigonometric KZ equations [^S], ^ and elliptic KZ-Bernard equations ||2^, for the 
six-vertex model and eight- vertex model, respectively. 
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